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a.  ru*rost  or  arm  ram 

To  compute  the  drawdown  time  (hrs)  of  a  reservoir  from  a  given  upper  pool 
elevation  to  a  specified  lower  pool  elevation  given  either  a  constant  inflow 
(may  be  zero)  or  an  inflow  hydrograph 
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•  .  PROGRAM  SPECIFICATIONS 

SEE  FOLLOWING  PAGES. 


C.  METHODS 

The  program  is  written  in  G635  time-share  series,  FORTRAN  IV,  and  is  part  of  a 
Conversationally  Oriented  Real-Time  Program-Generating  System  (CORPS).  The 
program  consists  of  a  main  program  and  a  subroutine.  The  main  program  handles 
all  I/O  requirements  and  computations.  The  subroutine  performs  all  writing 
and  file  requirements  to  save  the  output  for  graphics  and/or  other  use  if 
requested. 


D.  EQUIPMENT  err  AILS 

The  program  was  developed  and  is  operational  on  the  WES  G635>  Vicksburg,  MS. 
It  is  also  operational  on  HIS  66/80,  Macon,  GA,  and  Boeing  CDC,  Seattle,  WA. 


I.  ikftiV-MTNJf 

Input  requirements  are:  Either  constant  (may  be  zero)  or  inflow  hydrograph; 
time  in  hours,  flow  rate  in  acre- ft/hr,  loss  coefficient  of  outlet  conduit, 
initial  and  final  gage  heights  of  reservoir  water  surface  in  ft  above  outlet 
intake  invert,  and  coefficients  of  power  curve  to  data  of  gage  height  (ft) 
versus  reservoir  storage  (acre-ft).  Output  includes  given  data  and  the 
values  of  inflow,  outflow,  storage,  and  gage  height  of  water  surface  for 
each  time  step. _ _ _ _ _ _ _ _ 


r.  ADDITIONAL  REMARKS 

Complete  documentation  of  this  program  is  available  from  the  Engineer 
Computer  Programs  Library,  Technical  Information  Center,  WES. 
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B.  PROGRAM  SPECIFICATIONS: 

LANGUAGE:  ANSI  FORTRAN  (FORTRAN  IV ) 

Solution  Requirements:  The  run  command 

RUN  WESLIB/CORPS/H2250.R 


and  the  inputs  defined  in  (E). 

Method  of  Analysis:  Numerical  solution  of  the  level  pool  routing- 
lumped-form  of  the  continuity  equation. 

Core  Requirements  G635:  13  K  words 

External  Storage: 

1.  If  input  is  to  be  entered  for  an  inflow  hydrograph  and  the  user 
has  no  existing  hydrograph  file,  the  program  will  create  and  save 
the  user  named  hydrograph  file.  The  file  limits  will  be  from  6k0 
to  12K  G635  words. 

2.  If  the  output  is  to  be  saved  for  graphics  and/or  other  use,  then 
the  program  will  create  and  save  the  user  named  output  file  with 
limits  from  6kf6  to  12  K  G635  words. 

Restrictions :  None 

General  Equation:  Continuity  equation  discretized  over  a  time 
interval  (At) 
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*  y  -  2(oi 


where  I  ,  0  ,  and  V  are  inflow,  outflow,  and  volume,  respec¬ 
tively,  and  the  subscripts  1  and  2  indicate  the  start  and  end 
of  the  time  interval  (At). 
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REF :  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  722-F3-R0-2DO 

2.  TITLE:  H2250  -  Reservoir  Storage  Drawdown  Time 
3-  REVISION  LOG:  N/A 

U.  PURPOSE  OF  PROGRAM:  To  compute  the  time  required  to  draw  down  a 
reservoir  from  some  upper  pool  elevation  to  a  lower  pool  elevation  given 
either  a  constant  inflow  (may  be  zero)  or  an  inflow  hydrograph. 

5.  STEP  SOLUTION: 

a.  Input .  The  input  variables  that  must  be  supplied  upon  request 


are: 


(1)  Type  of  inflow.  Either  constant  (may  be  zero)  or  an  inflow 
hydrograph.  If  constant,  the  value  must  be  supplied  (acre- 
ft/hr),  however,  if  a  hydrograph  is  used,  a  data  file  must 
either  be  created  or  a  data  file  from  a  previous  run  may  be 
used.  If  a  data  file  does  not  exist,  the  program  will  fur¬ 
nish  instruction  on  how  the  user  may  create  the  file.  These 
instructions  include:  requests  for  data  file  name,  number 
of  points  on  hydrograph,  and  the  values  of  time  versus  dis¬ 
charge.  If  a  data  file  already  exists,  the  name  of  the  file 
is  requested. 

(2)  The  loss  coefficient  of  the  outlet  conduit  defined  by: 


where  K  is  the  loss  coefficient,  Q  is  the  discharge, 

A  is  the  cross-sectional  area  of  the  outlet  conduit,  and 
h  is  the  head  or  gage  height  (ft  above  intake  invert),  and 
g  is  the  acceleration  of  gravity.  See  program  H2251, 
which  may  be  used  to  compute  the  loss  coefficient  of  an 
outlet  conduit  from  prototype  measurements  of  drawdown. 
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(3)  The  cross-sectional  area  of  the  outlet  conduit,  ft  . 

(It)  The  initial  gage  height  of  the  water  surface  at  which  draw¬ 
down  begins  (feet  above  datum  used  to  determine  loss  coef¬ 
ficient  in  (2)  above). 

(5)  The  final  gage  height  of  the  water  surface,  ft. 

(6)  The  coefficients  of  a  power  curve  fit  to  the  plot  of  gage 
height  (ft)  versus  reservoir  volume  (acre-ft).  A  straight 
line  is  fitted  to  the  log-log  plot  and  the  equation  has  the 
form: 

Y  =  CX15  (2) 


where  Y  is  the  reservoir  volume  (acre-ft),  X  is  the 
gage  height  of  water  surface  (same  datum  as  in  (2)  and  (M 
above)  and  C  and  D  are  the  coefficients  required  by  the 
program.  If  one  straight  line  does  not  fit  the  data  over 
the  range  of  gage  heights,  straight  line  segments  of  best 
fit  should  be  used.  If  the  drawdown  limits  overlap  seg¬ 
ments,  the  problem  must  be  solved  in  steps,  where  the  ini¬ 
tial  and  final  gage  heights  for  each  step  lie  on  one  seg¬ 
ment  and  the  respective  C  and  D  coefficients  are  used 
for  that  step.  Thus,  the  final  gage  height  for  one  step 
would  become  the  initial  gage  height  for  the  succeeding 
step. 

(T)  Desired  type  of  output.  The  user  has  an  option  on  the 
amount  of  printout.  Either  a  detailed  printout  of  gage 
height,  inflow,  storage,  and  outflow  for  each  time  step  or 
a  printout  of  these  parameters  for  the  initial  and  final 
time  steps  are  available.  Furthermore,  for  the  detailed 
printout,  the  user  may  specify  the  time  step  increments. 

b.  Theoretical  Formulation:  The  level  pool  routing  technique  is 
used  to  describe  the  movement  of  water  through  the  reservoir.  The  pro¬ 


cedure  follows  from  the  lumped  form  of  the  water  continuity  equation. 


I  - 


(3) 


where  I  is  the  inflow,  0  is  the  outflow  and  V  is  the  reservoir 


storage  or  volume. 
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Equation  3,  discretized  over  a  time  interval  At  is 


f  *  V  '  f  <°1  +  °2>  * 


V  -  V 

_2 _ J. 

At 


(M 


where  the  subscripts  1  and  2  indicate  the  start  and  end  of  the  time 
interval,  respectively.  Equation  h  can  be  expressed  as 


Define 


+  V 


N 


V_ 

At 


+ 


0. 
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Equation  5  reduces  to 


N2=Nl-°l+2  (I1+V 


(5) 

(6) 


(T) 


The  data  consists  of  the  inflow  hydrograph  and  of  a  relation  between 
outflow  and  storage.  The  inflow  hydrograph  is  discretized  in  time  to 
yield  a  series  of  I  values. 

c .  Steps  in  Program  Execution. 

(1)  Input  variables  defined. 

(2)  Given  the  initial  gage  height  (head  on  outlet  conduit)  the 
initial  value  of  outflow  is  computed  by 


Q  =  0.0826UA' 


(8) 


where  Q  is  the  outflow  in  acre-ft/hr,  A  is  the  cross- 
sectional  area  of  the  outlet  conduit  (ft^),  g  is  the 
acceleration  of  gravity  (32.2  ft/sec^),  hq  is  the  initial 
gage  height  (datum  is  taken  at  the  outlet  intake  invert) 
and  K  is  the  outlet  conduit's  loss  coefficient. 


(3)  The  initial  storage  or  volume  of  the  reservoir  is  computed 
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V  =  Ch£  (9) 

where  V  is  the  storage  in  acre-ft,  h  is  the  initial 
gage  height,  and  C  and  D  are  the  coefficients  of  the 
power  fit  curve  of  gage  height  versus  storage. 

(4)  The  initial  N  value  (N^)  is  computed  with  equation  6. 

The  time  step  (At)  used  in  the  computations  is  set  at 

1  hour  for  the  conditions  of  constant  inflow  and  it  is  set 
at  the  time  interval  of  the  inflow  hydrograph  for  the 
hydrograph  condition. 

(5)  Given  the  inflow  hydrograph  (I  values),  the  updated  N 
value  (Ng)  is  computed  with  equation  7 • 

(6)  With  N2  known,  solve  for  h2  by  noting  that 

ch2  rrr  1/2 

"a'lr  +  °-0826AV/(h  2)  (10) 

or 

ah2  +  bh2/2  -  N2  =  0  (11) 

where  a  =  0.0826A*/^  and  b  =  ~  . 

\  K  At 

Equation  11  is  solved  by  a  first  order  Newton  technique, 

i  .  e .  , 


F(h)  =  ah°  +  bh2/2  -  N2 
and  the  derivative  with  respect  to  h  , 


(12) 


F'(h)  =  aDh^'1  +  |  bh‘1/2  (13) 

Initially  the  value  of  hg  is  assigned  the  value  of  h^ 
and  by  iteration,  i.e., 
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h  -  -F^}- 
h2  F’(h) 


(no 


u 


H2250 


the  new  approximation  (h^)  is  obtained.  Iteration  using 
equations  12,  13,  and  1^  is  continued  until 

<  0.01  (15) 


or  for  200  iteration. 

(7)  With  hg  determined,  the  values  of  storage  (V2)  and  out¬ 
flow  (O2)  are  computed  by  equations  9  and  8,  respectively. 
Set  V  =  and  0^  =  0^  and  repeat  for  the  next  h  . 

(8)  Termination  criteria  is  when  the  value  of  h?  is  less  than 
the  lower  limit  of  gage  height  supplied  as  input  or  when 
the  outflow  equals  the  inflow. 
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PART  II:  COMPUTER  FUNCTIONAL  DESCRIPTION 

1.  REVISION  LOG:  N/A 

2.  FUNCTIONAL  FLOW  CHART: 
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_  _  _  _ entry  point  into  H2250G  to  re- 
j  set  to  original  values  and  detach 

| H22501  |  | 

- _ _ 

reset  to  original  inter¬ 
val  and  output  type 
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3.  EQUIPMENT  AMD  OPERATING  SYSTEM:  The  program  was  developed  on  a 
0635  time-share  system  in  which  input/output  equipment  consisted  of  a 
Model  33  remote  teletype.  It  is  now  operational  on  the  WES/G635,  Vicks¬ 
burg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 

^ .  INPUT  REQUIREMENTS:  All  required  inputs,  except  for  the  case  where 
the  user  has  an  existing  hydrograph  input  data  file,  are  entered  via  the 
user's  time-share  terminal  device  in  free-field  format. 

5-  SECONDARY  STORAGE  INPUT:  The  following  is  the  format  setup  for  the 
user's  hydrograph  input  data  file.  If  the  user  has  no  existing  hydro¬ 
graph  file,  then  the  program  will  create  and  save  the  user  named  hydro¬ 
graph  file.  All  line  numbers  are  three  digit. 

Line  1  FORMAT (LX,A5, IX, 13) 

Equal  to  H2050  followed  by  the  number  of  points  (NP)  in  the  hydrograph. 
Must  always  have  H2050  since  it  is  used  by  WESLIB  routine  TACHFILE  to 
verify  that  the  file  is  an  input  data  file  for  program  H2050 
100  H2050  19 

Line  2  FORMAT(3X,2(lX,F10.2) ) 

This  format  controls  the  remainder  of  the  data  file.  Each  line  will  con¬ 
tain  one  value  of  time  and  one  value  of  discharge.  The  number  of  lines 
will  be  equal  to  the  number  of  points  (NP)  as  given  for  the  hydrograph. 

102  0.00  0.00 

10L  1.00  6.6l 


2NP+100 


18.00 


0.00 
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6.  INPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the  input 
variables  in  program  H2050. 

2 

AC  -  outlet  conduit  cross-sectional  area,  ft 

CD  -  loss  coefficient  of  outlet  conduit 

Cl  -  constant  inflow,  af/hr 

COEFFA  -  coefficient  of  power  equation  for  gage  height  (ft)  versus 
reservoir  storage,  af/ft 

COEFFB  -  coefficient  of  power  equation  for  gage  height  (ft)  versus 
reservoir  storage,  af/hr 

FILEK(l)  -  8  character,  dimension  max  2;  name  of  input  data  file 

HI  -  initial  gage  height  of  reservoir  water  surface,  ft 

HF  -  final  gage  height  of  reservoir  water  surface,  ft 

IPR  -  time  step  increment  for  printout 

NP  -  number  of  points  on  inflow  hydrograph,  <100 

NT  -  type  of  inflow,  0  -  constant;  1  -  hydrograph 

OUTPUT  -  integer;  type  of  output,  1  -  detail;  2  -  final  conditions 
only 

Q(l )  -  dimension  max  100;  inflow  values  of  hydrograph,  af/hr 

T(l)  -  dimension  max  100;  time  values  of  hydrograph,  hr 

7.  OUTPUT  DATA  DESCRIPTION :  The  following  names  are  used  for  the  out¬ 
put  variables  in  program  H2050. 

H2  -  gage  height  of  reservoir  water  surface  at  end  of  time  step, 
ft 

HI  -  initial  gage  height  of  reservoir  water  surface,  ft 
01  -  inflow  at  initial  time  step,  af/ft 

02  -  inflow  at  end  of  time  step,  af/hr 
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Q2  -  outflow  at  end  of  time  step,  af/hr 

QI(l)  -  dimension  max  100;  outflow  at  initial  time  step,  af/ft 

TT  -  time  at  time  step  n  for  n  =  initial  time  to  drawdown 

time 

V0L1  -  reservoir  storage  at  initial  time  step,  af 

V0L2  -  reservoir  storage  at  end  of  time  step,  af 

8.  PROGRAM  MESSAGES:  Return  from  program  message  is  to  RERUN  (another 
set  of  conditions)  or  STOP. 

NEWTON  METHOD  DID  NOT  CONVERGE  IN  200  ITERATIONS  -  CHECK  INPUT  DATA 

9.  VARIABLE  DEFINITIONS; 

AA  -  temporary  variable  used  in  computing  the  gage  height  of 

the  reservoir  water  surface  at  end  of  time  step 

2 

AC  -  outlet  conduit  cross-sectional  area,  ft 

AN1  -  N  value  at  beginning  of  time  step  (N^) 

AN2  -  N  value  at  end  of  time  step  (N^) 

B  -  temporary  variable  used  in  computing  outflow 

BB  -  temporary  variable  used  in  computing  the  gage  height  of 

the  reservoir  water  surface  at  end  of  time  step 

CD  -  loss  coefficient  of  outlet  conduit 

Cl  -  constant  inflow,  af/ft 

COEFFA  -  coefficient  of  power  equation  for  gage  height  (ft)  versus 

reservoir  storage,  af 

COEFFB  -  coefficient  of  power  equation  for  gage  height  (ft)  versus 
reservoir  storage,  af 

DELT  -  time  step,  hr 
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DH  -  difference  between  previous  value  of  gage  height  at  end 

of  time  step  and  the  value  of  the  new  approximation, 
equation  15,  ft 

DQ  -  difference  between  outflow  and  inflow,  af/hr 

FH  -  functional  relation  (equation  12)  in  terms  of  gage  height 

at  end  of  time  step 

FHP  -  first  derivative  of  FH,  equation  13 

FILEK(l)  -  8  character,  dimension  max  2;  1=1,  name  of  input 

data  file;  1=2,  name  of  output  data  file  for  graphics 
and/or  other  use 

2 

G  -  acceleration  of  gravity,  32.2  ft/sec 

HI  -  new  approximation  of  gage  height  of  reservoir  water 

surface  at  end  of  time  step,  equation  lU,  ft 

H2  -  gage  height  of  reservoir  water  surface  at  end  of  time 

step,  ft 

HF  -  final  gage  height  of  reservoir  water  surface,  ft 

HFILE  -  5  character  name  of  program  (H2050);  passed  to  WESLIB 

routine  HACCT  for  bookkeeping;  passed  to  WESLIB  routine 
TACHFILE  to  verify  that  I/O  files  to  be  used  are  files 
for  H2050 

HI  -  initial  gage  height  of  reservoir  water  surface,  ft 

IPR  -  time  step  increment  for  printout ,  hr 

JKL  -  directs  return  from  WESLIB  routine  RERUN  to  desired  input 

read 

KFILE  -  file  code;  KFILE  =  1  ,  file  is  an  input  data  file;  =  2, 
file  is  an  output  file 

LQX  -  equal  1,  print  instructions  from  WESLIB  routine  RERUN; 

equal  3,  no  print 

LQ7.  -  equal  1,  execute  all  manual  input  cues  and  reads;  equal  2, 

call  WESLIB  routine  RERUN  and  enter  only  desired  inputs 

NP  -  number  of  points  on  inflow  hydrograph 
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NT  -  type  of  inflow;  0  -  constant;  1  -  hydrograph 

01  -  inflow  at  beginning  of  time  step,  af/hr 

02  -  inflow  at  end  of  time  step,  af/hr 

OUTPUT  -  integer;  type  of  output,  1  -  detail;  2  -  final  conditions 
only 

Q1  -  outflow  at  beginning  of  time  step,  af/hr 

Q2  -  outflow  at  end  of  time  step,  af/hr 

Q(I)  -  dimension  max  100;  inflow  values  of  hydrograph,  af/hr 

T(l)  -  dimension  max  100;  time  values  of  hydrograph 

TT  -  time  at  time  step  n  for  n  =  initial  time  to  drawdown 

time 

V0L1  -  reservoir  storage  at  beginning  of  time  step,  af 

V0L2  -  reservoir  storage  at  end  of  time  step,  af 

ZZZZZ  -  2  character;  equal  RE,  rerun;  equal  ST,  stop 

10.  EXAMPLE  CASE:  User  has  an  existing  hydrograph  data  file. 

Input:  Type  of  inflow  (NT)  =  1 

File  name  (FILEK(l))  =  HYDRO 

Loss  coefficient  outlet  (CD)  =  16.2 

Cross-sectional  area  (AC)  =  16  ft2 

Initial  gage  height  (Hi)  =  60  ft 

Final  gage  height  (HF)  =  10  ft 

Coefficients  of  power  curve  (COEFFA)  =  .0^35; 

( COEFFB )  =  2.57^ 

Type  of  output  (OUTPUT)  =  1 
Time  interval  (IPR)  =  12  hr 

Save  output  for  graphics  and/or  other  use  =  NO 
List  of  HYDRO 

102  m22S0  19  120  9>f*3  14.38 

e,  0.  122  *0.«®  13.22 

10a  1.00  6.61  *2*  11.00  11.57 

106  2.00  la. 05  126  12.00  9.92 

103  3.00  21. «9  I23  13.00  3.  26 

no  a. 03  2a.  79  *30  la. 00  6.61 

112  5.00  22.  31  132  15.00  a.96 

Iia  6.00  19.83  13a  16.00  2.40 

116  7.00  18.  13  136  17.00  0.33 

113  3.33  16.  53  138  13,00  0. 


H2250 


Output : 


H2250  -  pESEPUOI P  5  TO  PAGE  EUACUPATIOM  ”IME 

AA  -  ”VPE  OP  INFLOW*  EM”E®  0  TO®  COM«”AMT  INFLOVCMAY  BE  ?EPO>  O- 
PITE®  1  TOP  INFLOW  HvDPOGPAPH(MA*'  BE  FOLLOWED  WJTH  CONSTANT  INFLOW). 

-1 


CO  vO’»  HA"E  AM  INPUT  DATA  FILE  FOP  INFLOW  HYDPOGFAPH*  Y  OP  N! 

IF  MO#  »I  WILL  HELP  YO”  CPEATE  FILE. 

— v 

AF  -  FILE  NAME 
•HYDFO 

AG  -  LOSS  COEFFICIENT”  FOP  0U”LE”  CONDUIT  DEFINED  Pv> 

K-2.-G*H/(P/A)p*2*  vhEPE  K-LOSS  COEFF,  G“ACCELEFATtON  OF 

GPA"I  ’V  ,  P*  DI  « CH  APGE  AMD  A-C®0  FF-  <EC”1  OMAL  APEA  OF  0"TLET  CONDUIT. 

-16.2 

AH  -  CPO'S- SECTIONAL  APEA(FT**2>  OF  OUTLET  CONDUIT. 

•16 

AI  -  INITIAL  GAGE  HEIGHT  AT  WHICH  DPAVDOWM  BEGIMFt  FT.  ABOVE 
DATUM  USED  TO  DETEFMINE  OU”LET  COMD”lT  LOSS  COEFF). 

-60 

AJ  -  FIMIAL  GAGE  HEIGHT  TO  VHI CH  PESEPVOIP  IF  TO  BE  DPAVN  DOWN 
CFT.  ABOVE  DATUM  VFED  TO  DETEPMIME  OV”LET  COND"IT  LOFS  COEFF). 

-10 

AX  *  COEFFICIENTS  OF  "0”E®  C"»’’E  FIT  TO  GAGE  HEIGH”(  FT)  VS  PESEP*TOIP 
VOLUME!  AC—E-FT)  •  FOPM  OF  EIV'ATION*  Y-A*X*«P*  WHEPE  v  EOUAL c 
VOLT,ME,  Y  EO”AL?  GAGE  HEIGH”  AMD  A  AMD  B  APE  THE  COEFFICIENTS. 

E?tep  a  -HEM  B  SE®E®ATED  WITH  A  COMMA. 

-.0435, 2. 574 

AL  -  DESIPED  type  OF  9"*"''*  AS  FOLLOWS* 

DITEP  1  -  DETAILED  P PINTO””  OF  GAGE  H”.,  INFLOW,  STOPAGE 
AMD  OUTFLOW  FO*»  EACH  IMCPEMENTED  ”IME. 

EM”EP  2  -  LIMITS  0VT?”T  TO  INITIAL  AMD  FINAL  CONDITIONS  ONLY. 

-I 

Aft  -  TIME  IMTSPUAL  FOP  PPIMT0”T. 

MOTE*  MUF”  PE  E'V’AL  ”0  MULTIPLES  OF  I  FOP  CONSTAM”  INFLOW  OP  MVL”I?LE? 
OF  time  INTEPVAL  OF  THE  INFLOW  HYDPOGPAPH. 

AN  -  *AVE  OUTPUT  FOP  GPAPHICS  AND/OP  OTHEP  ”SE 
•MO 


TIME 

GAGE  H”. 

INFLOW 

STOPAGE 

OUTFLOW 

(HP) 

<  FT) 

<  AF/HP) 

<  AF) 

(AF/HP) 

0. 

60.00 

0. 

1642. 29 

20.42 

12.00 

59.20 

9.92 

153  6.  50 

20.28 

24.00 

55.99 

0. 

1374. 19 

19.73 

36.00 

52.09 

0. 

1141.59 

19.03 

40.00 

47.36 

0. 

917.39 

18.24 

lh 


-5' 
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60.00 

43.  13 

0. 

704. 37 

17.  32 

72.00 

37.33 

0. 

502.34 

16.  23 

34.00 

31.63 

0. 

316. 12 

14.83 

96.00 

23.  64 

0. 

149. 34 

12.82 

103.00 

10.07 

0. 

1  6.  62 

8.  37 

109.00 

7.85 

0. 

8.75 

7.39 

t»rpf<«I  OP 

■  CT^p 

STO° 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineer  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  III:  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE:  H2250  -  Reservoir  Storage  Drawdown  Time 

3.  PROGRAM  SOURCE  LISTINGS:  See  pages  17-22 

U.  NUMERICAL  AND  LOGICAL  ANALYSIS: 

First  order  (Newton)  fixed  point  iteration  technique  is  used  to 
determine  gage  height  of  reservoir  water  surface  at  end  of  time 
step 

5 .  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  None 

6.  MISCELLANEOUS :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational  on  the  WES 
G635,  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 
The  source  listing  on  page  17  contains  the  first  line  run  command  and 
brief  for  H2250.  This  first  line  run  command  runs  the  binary  H2250B  of 
the  source  listing  on  pages  18-22  (FORTRAN  source  of  H2250)  and  attaches 
the  WESLIB  routines  HACCT,  TACHFILE,  and  RERUN. 
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H2  2C;0 


H2250 


02/06/79 


0001*#PPM  VE*LI  P/  CO  »'t'p/H 22503.,  pi  VEFL  I  r/PEPPM,  pJ  VE^LI  9/~ACHF!  LE,.  P* 

0032*  0'iE'Ll  P/HACC  T*  p 

38  0  0  6  2  THI?  PPOGPAH  COMPETE®  THE  DPAVDO'FJ  TIME  (HP)  FOP  A  PEcEP"OIP 
03  1  0  64vj  ~~r{  CO‘J  eTAM T  I MFLO"  (  AF/H°)  OP  IMFLO”  HYDP9GPA',H  FnECI  FI  ED  VI 'H 
38  2  0  6  2-HE  GI  ’’EM  GAGE  HEIGHT  (F_)  ”EFF"F  '"“OPAGE  PELA”rIO*l  ( AF)  .  OTHEP 
0330  59im®'Tpf  PEG"!  PED  APE  THE  C®0 FF- FEC~: OM  AnEA  (FT**2>  AMD  LOFP 
0340  60COEFFI  Cl  EM-  OF  —HE  0"TLET  COMP’HT  AMD  TH  E  L1M1  (Ft>  OF  THE 
38 50  24"ATP  pv^FAGE  EDE,,A-rI  OM. 

036C  63  THE  0”Tp”t  COMF1  FT  OF  A  ’'APLE  OF  TIME  CH®>  ’'EPF"®  GAGE  HEIGHT 
0870  65<  PT)  ,  I  “IFLO”  (  A  F/H  F)  *  ®T0  PAGE  (AF)  AMT  0*’-FL0«  (AF/HP)  A  F  '“EC! FI  ED 
0380  A75V  THE  "«£?'?  TIME  IMCFEMEMT  (HP)  FOP  PPIMTO"'-. 

0999*06f!MI FH 
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32/06/79 


H2250e 


0000  t*#p»iM  *»)  /  CO  P°S/H2250BC  VOGO  ) 

10000  DIMSTIFIOM  f(  100),  oi  C  100) 

10010  COMMOM  /C77331/KFILE,  FJLEXC 2), /MAIV/LCT.LeM 
10020  CMA^AC-E®  F!LEK*8 ,  H FI L E*  5,  * ZtZZ*  2 
10030  IMTEGEP 

10040  L«T*  l>LF'*  1'  MFI  LE*5hh2250)*ICO>1*  W V~C“ I 
10350  "PIMT  13333 

10060  10033  FOPMATC /"H2250  -  »E'EP"0!P  FTO PAGE  E''ACT'PAT!OM  TIME"//) 

10070  10034  CALL  HACCTCMF1LE) 

10080  GO  TOC  10040,  10346),L(JJ 
10090  10040  *>PIM*  I 

10100  1  FOPMATC  "AA  -  T^PE  OF  IMFLOV*'  ENTEP  0  FOP  COM'*AVT  1  MFLO'/c Ma ■*  PE 
10110«TE"O)  OP’V’EVTEP  1  FO r  INFLOW  HV DPOGPAPH C MAT  PE  FOLLOVED  VITO  COS' 
101204TAMT  IMFLOV)  .  "> 

10130  10030  »EAD,MT 
10140  IFCMT.ME.  3)  GO  TO  2 
10150  I FCMCOM. EO.  1 )  GO  TO  10100 
10160  GO  "OC 10100, 10346), LOT 
10170  10100  ppjriT  A 

10130  4  'OPMATC "At?  -  CONFTAMT  IMFLOV  VALVE  IV  ACPE- FT/HP.  IF  MO  IMFLOV 
101904-e*ITEP  0.  0. 

10200  10130  PEAD, Cl 
10210  MC0M“2 
10220  CI< l)*Ct 
10230  BELT* 1.0 
10240  GO  TO  30 

10250  2  GO  TOC  13175,  10205), Lr’? 

10260  10175  ipiMT  6 

10270  6  «-0®‘!ATC /"DO  v0"  >IA"E  A’l  IV'»,,T  DA-A  FILE  FOP  IMFLOV  HvPo0G'!A®m7 
102304  -  OF  MfVlF  *J0,  VE  VILL  HELP  ""V  CFEATE  FILE.") 

10290  CALL  AM'VEFC J5CP, $1 50) 

10300  10205  I  FCMH'*.  EG.  1  )  L0E*1 
10310  MTT'-S 

10320  GO  ■’DC  1  50,  10846), LOT 
10330  150  PPINT35LOT*  l 

10340  3  FOPMATC "AC  -  ME V  IMP"T  DATA  FILE  MAME.") 

10350  10230  r>EAO  9,  FILEKC  1 )  «LOE*  l 
10360  9  F0«»MATCA8) 

103  10  GO  TOC  10250,  10346), LO”. 

,i33e  10250  vopjT  10 

10390  10  FOPMATC  "AD  -  MvMPEP  OF  POIMTF  FOP  HY DPO G PAPH .  " ) 

10400  10270  PEAD,M« 

10410  GO  TOC  10230,  10346), L*Z 
10420  10230  r«?IMT  11 

10430  II  FOPMATC  "AE  -  PAI PJ  OF  TIMECH®?)  «»S  OC  ACPE-FT/HP)  OF  THE  Hvd®OGP 
104404ArH) "/" 'EDEFATE  ALL  ”AL»»E'  VJ “H  COMMA?.  MAYIM”M  MO. “100  AMD  TIME"/ 
104504"  I  M7E?'’AL  MV  FT  PE  COM?TAMT.  FOLLOVJMG  HY  DPO  G  PAPH  "HE  PPOG»AM  A??VM 
104604E'  A"/"COM 'TA.VT  IMFLOV  EfH’AL  TO  LA5T  C.  VAL’’E  OF  THE  HYPDO G®APH . ") 
10470  10330  PEAD, CTCI)»CICI),I*1, VP) 

10480  KFILE-1 
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M2250'-  COMT. 


10490  CALL  TACHFILFCHFILE.  J  l  1  734) 

10500  ””I -ECKFILE. 270)  M", ( I  *  2+ 1 00,  t< I ) » PI ( I) *  I “ 1 ,MP) 

10513  CALL  DETACH C  KFILE,  »  > 

10520  27 0  r7'"lA'-(  "100  E"ACU  *M  3/ f  X  3,  2(  ly,  Fl  0.  2)  )  > 

10530  CD  ~0  20 

10540  503  !  EC!"".  ME.  1 )  CO  -r0  20 
10550  "»!»!’, " Ar  -  FILE  MAME" 

10560  13410  or/\p  9#FILEII<1> 

10570  :<FILE=1>‘LCT  =  2 

10530  CALL  '■ACHFILE(HFILE>  II  l  734) 

10590  c  E v  I M  D  KFILE 
10600  IECLOr.EA.l)  LnT=l 

10610  PEADCKF1LE,  540,  EMD=54S)  MP,  C  T(  I  )  ,  Cl  C  I  >  t  I  =  1  ,  MP) 

10620  540  FOFMATC  10-'<  !  3/c  3Xj  2<  lv,  pl0.  2)  >  > 

10630  545  COMTIMVE 

10640  CALL  DETACH C KFILE*  *  ) 

10650  20  DELT=  TC  2)  -  TC  1 ) 

10660  30  I FCLOV. En.  1 )  GO  TO  35 
10670  L0T=2 

10630  35  GO  •'Ot  leSIS,  10346), LOT 
10690  10515  pPIH'-  13 

10700  13  FOPMATC  "AG  -  LOFF  COEFFI  Cl  EPJ  T  FO"  0"n.E*  COMtH'I  T  DEFIMED  FT!"/ 
10710i"K=2.*G*H/CO/A)**2J  "HEPE . K=LO FF  CO EFF>  G= ACCELEPA’'! O'J  OF’VGPAVITT 
107204, n=DI eCHA°GE  AMD  A=  CPO FF- 5ECTI OMAL  APEA  OF  0"TLET  COMD"IT."> 

10730  10550  pEAD, CD 

10740  GO  TOC  12562,  13346), LOT 

10750  10560  opIJJT  14 

10760  14  FO  "MATC  "AH  -  CFO FF- eECTI OMAL  APEACFT**2)  OF  0,'TLE*'  COMD'MT.") 

10770  10580  pEADjAC 

10780  GO  TOC  10590,  10846), LOT 

10790  10590  nplMT  16 

10800  16  FO  "MATC  "AI  -  IMITIAL  GAGE  HEIGHT  AT  WHICH  DPAVPOVM  PEGIMFCFT. 

108104  AS0”E"/"DATT71  ”FED  T0  DETE"MIME  OUTLET  C0MD1JIT  LOP?  COEFF..") 

10320  10620  pEAD.HI 

10830  GO  tqc 10630, 10346), LOT 

10840  10630  PPIMT  17 

10350  17  FO  PMATC  "AJ  -  FIMIAL  GAGE  HEIGHT  TO  VH  I CH  PEFEP''OIP  IF  TO  BE  DPA 
103604VTJ  COVM"/”C  FT.  ABO’'E  Datum  tjceD  TO  DETEPMIME  0”TLET  COMD"IT  LOFS  C 
103  7040  EFF,  ") 

10330  1B660  PEAD.HF 

10390  GO  -0C 10670, 10846), LOT 

10900  10670  ppiMT  13 

10910  18  FOPMATC  "AK  -  CO  EFFI C I  E*IT'!  OF  POVEP  CPP"E  FIT  TO  GAGE  HEIGHTCFT) 
109204  ”F  PEFEP”0 1  P"/'"'0L’P1EC  ACPE-FT) .  FOPM  OF  EP”AtiOM!  v“A*,C**P'  VHEPE 
109304  ”  E01'ALF"/"'TOL’F1E<  M  EP"ALF  GAGE  HEIGHT  AMD  A  AMD  B  APE  THE  COEFF 
1094041 C!  E")Tt:.  "/"EMTEP  A  THEM  P  'EPEPAtED  VI  TH  A  COMMA.”) 

10950  10720  PEAD* COEFFA,  COEFFP 
10960  GO  TOC  10730,  10346), LOT 
10970  10733  p«IMT  19 

10930  19  FOPMATC  "AL  -  DEFIPED  •p'PE  OF  OT!TP”T  AS  FOLLOWS!  ”/5k*  "EDJTEP 
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H2250e  CO'IT. 


109904  l  -  DETAILED  ?P1MT0”T  OF  CAGE  XT..  IMFLOV.  ?TO°AGE"/ 1  6*',  "AMD  0"TF 
1 1 0004LO V  TOP  EACH  IMCPEMEMTED  TIME.  "/5r>  “E'l^EP  2  -  LIMIT?  OMTP"T  TO  IMI 
1 1 0  1 04’’!  AL  AMD  PIMAL  COMDfMOMe  OMLV.  ") 

1 1020  KO’”-*  I 

11030  10780  pEAD. 0"TPVT 

11040  GO  TOC  107S7, ie346).L OT 

U05G  10787  I  *,(0,,TT>*,',".ME«  1 )  GO  TO  29 

11060  I  FCKO"T.  E'’.  1 )  GO  *0  10300 

11C70  GO  TOC  10300,  10846).LOT 

11030  10300  opIMT  12JK0,,t*2 

11090  12  FOPMA-rC”AM  -  TIME  IMTEF*>AL  FO*1  PrtMTO,"r.  ”/”M0  WS*  M"**,‘  BE  ECVAL 
111004  T)  M»’LTI’>LEE  OE  1  FO^  COM«TAMT  INFLOW  0®  MfLTI  PLEe"/.  6r.  “0 E 
111104  TIME  INTERVAL  OF  THE  INFLOW  MYDFOGFATTl. "> 

11120  108  40  ®EAD. I"® 

11  130  29  GO  TOC  10350,  10846), LOT 

11  140  10850  0*1  M-,  "All  -  SA"E  O'lTFT'T  FOP  GFAPHIC?  AMD/OP  OTHEP  *»SE" 

11150  CALL  AM?'?EC(  S  1 1  1  S3,  *  l  1  1  58) 

11160  11153  M®A"E»1 
11170  PPIMT. "FILE  MAME“ 

11  130  KFILE»2 

1119e  11156  PEAD  9, FILEKC  2) 

11200  GO  TO  I  1  160 

11210  11158  M  ?ATfE“0 

11220  11160  GO  TOC  11200,  108 46). LOT 

11230  10846  KHK*14 

11240  CALL  “EF'TICKKX.Lor.JKL) 

11250  GO  TOC  10080,  10130,  10230,  10270,  10330,  10410,  10550,  10580,  10620,  10660, 
11260410720,  10780,  108  40,  11156, 1 120C),JKL 
11270  1 1200  n®IMT  31 

11280  31  EOr'lATC//.  2^,  "TIME  GAGE  HT.  INFLOW  FTOPAGE  Ot'TFLOV"/" 
11 2904  C  HP)  ",  7t,  "C  FT)  ".  Sr,  "<  AF/H  P>  ”/  3r.  “C  AF>  ",  5r.‘  "C  AF/HF)  ”>  /) 

11300  MGF*0 
11310  11225  H 2“ HI 
11320  K*0 
11330  G* 32, 2 
11340  I«l 

11350  P»0. 08264* AC*  POPTC  2.  * G/CD) 

11360  TT»0.0 
11370  0 1"P* ?®PTCH! ) 

11380  'COL  1*C0EFFA*HI ** COEFFB 
11390  IFCMG'.EO.l)  GO  TO  11322 
11400  PPIMT  34, TT.HI, 01 C • ). "OLl.O • 

11410  GO  TO  90 

11420  34  F0FMATCF7.2,  5r,  F6.2,  4T.F6.2,  3v,F7.2,4r,F6.2) 

11430  11322  PPI TEC  2,  1  1 323)  TT.HI.OICl). POL  I  *  0 l 
11440  11323  F0P51ATC  3f3.2.  F10.2.  f8.2) 

11450  40  AMl“"0Ll/DELT*0l/2. 

11460  K”K*1 

11470  IFCMT.ME. 0)  GO  TO  45 
11480  Cl”CI 
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11490  n  2=  Cl 

11500  45  I  FCMt.*]F.  1 1  GO  TO  50 
11510  IECI.GE.H'')  GO  TO  51 
1 1  520  r>l-=r»i(i5 
11530  r2=C! C I + 1 ) 

11540  CO  TO  53 
11553  51  Pl  =  CI(M”5) 

1 1  560  ',2=r'J  CMP) 

11570  50  AM  2=  ATI  1  -0  1  +  C  0  1  +  02)  /  2. 

11533  AA=COEEEA/DELT 
11593  pr=P/2. 

11600  KH=0 

11610  55  EH  =  H  2*  *  CO  EEE3+  C  FB/AA)  *B0PTCH2)  “AM  2/ A  A 

11620  E”H  =  C0EEEB*H2**C  COEFEB-  1 . 0)+CBB*0.  5/AA)  /  SCP'I'CHE) 

11630  H 1  =  H 2- FH/ EnH 
11640  PH=AB'(Hl-H2) 

11650  KK=KK+ I 
11660  H  2=H 1 

11670  I  rC  KK.  GT*  200)  GO  TO  60 
11630  I  =C  PH.  GT.  c.  0  1 )  GO  TO  55 
11690  GO  TO  65 
11700  6C  p*3 1 M T  61 

11710  61  FOFMATC/"MEVTOM  ME'-'HOD  DID  MOT  COM”EPGE  IM  200  ITE^ATIOMS  - 
1 1  720^6 CHECK  IM”M-  PATA.") 

11730  GO  TO  999 
11740  65  TT=TT+DELT 
11750  02=P* crnTf h2) 

11760  ’'OL 2=  CO EEEA+H  2**  CO  EEEB 
11770  pc= APEC A2-0 2) 

11730  I  EC  0"T?t’T.>JE.  1 )  GO  TO  66 
11790  KK^K/IP^+I^P-K 

11300  IFCKKl.L'r.0.OP.KKl.GT.0)  GO  TO  68 
1131B  67  I FCMG”. EO. 1 )  GO  TO  11696 
11320  np’MT  34, -'■t,h2,  Q2,  V0L2..  02 
11330  GO  TO  66 

11340  1  1696  "PI  TEC  2,  1  1  323)  TT»  H2,  o2,  MOL 2,  0 2 
11350  66  I FCMGC. EC. 1 )  GO  T0  11703 

11360  I  ECO GE.  2.  AMD.H2.LT. HD  PPIMT  39,  TT, H2»  q2,  ’I0L2,  0 2 
11370  GO  TO  1  1  710 
11380  11703  COMT!M”E 

11390  11710  I  ECH2.  LT.  HE.  OF.  DP.LT.  .  0  1  >  GO  TO  999 

11900  63  I  ECH2.  L^.HE.  OF.  DO.LT.  .  0  !  )  GO  TO  67 

11910  1=1+1 

11920  01=02 

11930  "0t,l  =  ”9L2 

119  40  GO  ”0  40 

11950  999  L0/  =  2 

11960  IECMC’.E0.  1)  CALL  H2  250  1  c  OHTF'TT+  I  PP,  1 1  1  73  4) 

11970  I  FC M ^A”E.  En.  1 )  CALL  H 22 50GC HEIL Z,  "T. 0"‘rP”‘r/  I  ”P+  MG”#  1 1  l  78 4,  J 1  1  225) 
1 1930  1  1  734  r»PIM'r  1  1  785 
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H2250P  COMT. 


11990  11735  FO"MAT<  //"PEPIPJ  OF  STOP") 

12000  11430  PEAD,TTTTT 
12010  I  F<  ZZZZZ  •  E".  2:ifF)  OO  TO  10734 
12620  I  F(ttt7T.  ,  jn.  2Tl"T)  00  TO  11430 
12030  ”r,I *IT>  "AMO  M"ST  BE  PE0""!  OP  STOP" 

12040  GO  TO  11430 
12050  11430  STOW 
12060  END 

13000  Sl’PPO  VT1ME  H2250QCHFILE.-  TT,  0"TP'.'T*  1  "P,  MGP, *, *  ) 

13010  COMMON  /C77301/KFILE, FILEKC 2) /MAZM/LOT, LOM 
13020  CHA°ACTEr  FILEK*8,HFILE*5 
13030  IMTEGEP  Ol'TP"T 
13040  KFILE=2 

13050  CALL  TACHFILECHFILE, *13220) 

13060  ND«TTJJIV"TTJJl2»l 

13070  13060  I FCTD.LE. 200)  GO  TO  13080 

13030  m2=jj2*2JMD="JD/M2J  GQ  TO  13060 

13090  13030  MSl  =  0'.'TP,'Ti'ls2=Ipp 

13100  I  PP=TJ2J  0"TP"T=  1 

13110  f (  M”/M  2)  *’12 

13120  IFCNV.GT.0)  ?JV*  I 

13130  MD=MD+ 1+JJV 

13140  pEVIMD  2 

13150  VPI TEC  2, 13120)  <ND,I«1,5) 

13160  13120  FO PMAT< "10  H2250  05  05  EDGE"/" 1 1  < 05( /) * < 3f8. 2# Fl 0. 2,  f3 . 2) > " 
13170«/"12  ",  I  3,  4<  1y,  x  3)/"13",  5<  ’•  l")) 

13160  VPITE<2, 13150) 

13190  13150  F0PMAT<"14  CTIPVE  DESI  GMATIOfIS  FOP  H2250  APE«”/"15  1»TIME 
1320042=GAGE  HT  3=  INFLOW”/"  1  6  4-STOPAGE  5“0’’TFL0W”/"1 7  ’PUTS  FOP  AB 

132 1 04O"E  "API  AFLES  APE!'V"18  HP«1  FT“2  AF/HP"3,  5  AF*4") 

13220  ’JGP»1 
13230  PET"PM  2 

13240  EflTpr  H22501(3»tTPT'T,  IPP,*) 

13250  I  PP“M  $2J  OUTPUT-?!  S 1 1  CALL  DETACH{2,#> 

13260  13220  PET"DM  1 
13270  EMD 
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